NACA TN 4315 ;qqnT 



NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 4315 





AN ESTIMATE OF THE FLUCTUATING SURFACE PRESSURES 

ENCOUNTERED IN THE REENTRY OF A BALLISTIC MISSILE 

By Edmund E. Callaghan 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 



Washington 

July 1958 


' i 

■IT- 



TECH library kafb, nm 



TECH LIBRARY KAFB, NM 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTIC^ 


TECHNICAL NOTE 4315 


AN ESTIMATE OF THE FLUCTUATING SURFACE PRESSURES ENCOUNTERED 
IN THE REENTRY OF A BALLISTIC MISSILE 
By Edmund E. Callaghan 


SUMMARY 

Calculations are made of the magnitude of the surface sound pres- 
sure levels that will occur on blunt-nosed vehicles during reentry into 
the earth's atmosphere. The results presented cover a wide range of re- 
entrance velocities and reentrance angles into the atmosphere and sizes 
and weights of the vehicle. 

The analysis presented is based on current subsonic results, which 
can be extrapolated to high supersonic speeds for a blunt-nosed body 
since the values of local Mach number behind the shock but outside the 
boundary layer are quite low. The analysis is limited to the particular 
point on the body where the maximum local sound pressure level is ob- 
tained. A constant value of the ratio of the surface root-mean-square 
pressure to the local dynamic pressure, which is based on current sub- 

>7 

sonic data, is assumed (4.5X10“°). 

The results indicate that surface sound pressure levels of the order 
of 150 decibels or higher are likely to be encountered. The time of ex- 
posure to these levels is found to be of the order of 25 seconds. 




INTRODUCTION 

The problems associated with the reentry of a ballistic missile are 
numerous. In some cases even the order of magnitude of the problem is 
unknown, such as the internal and external noise levels to which the nose 
cone is subjected. High external fluctuating pressures in the boundary 
layer may damage the external structure of the vehicle if the times in- 
volved are of long enough duration. On the other hand, high external 
sound pressure levels may result in high internal sound pressure levels 
(inside the cone) and may interfere with the proper action of internal 
mechanical or electronic devices. 
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The analysis presented herein was made to determine the maximum 
value of sound pressure level that would occur during the reentry of a 
hlunt-nose cone. 

Although there are very little data on the fluctuating pressures in 
a boundary layer even at subsonic speeds (refs. 1 and 2), the necessity 
of making an order-of -magnitude estimate of the sound pressure levels 
inside and outside the nose cone is obvious. The results reported herein 
represent an extrapolation of current knowledge and data in the subsonic 
speed range to the very high hypersonic speeds encountered by ballistic 
missiles. 

Such an extrapolation would appear to be possible for blunt-nosed 
vehicles, which are of current interest, since the Mach numbers around 
the body (behind the shock wave) are only slightly supersonic. Further- 
more, the measurement and use of acoustic data on the decibel scale (6 
db is a ratio of 2 in rms pressures) to a probable accuracy of no better 
than ±2 decibels is usually sufficient and, hence, errors involved in 
data extrapolation may not be any larger than the usual accuracy to which 
the final result is desired. 


SYMBOLS 

A area, sq ft 

a speed of sound, ft/sec 

Cjj drag coefficient 

e base of Naperian logarithims, 2.718 . . . 

g acceleration due to gravitational force, 32.2 ft/sec 2 

M Mach number 

m mass, slug/cu ft 

P total pressure, lb/sq ft 

p pressure, lb/sq ft 

P root -me an- square pressure at surface 

q dynamic pressure, (r/2)M 2 p 

R gas constant 
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Re Reynolds number 

SEL sound pressure level, db(re 2x10“ 4 dynes/sq can.) 

T temperature, °R 

t time, sec 

V velocity, ft/sec 

V E velocity at entrance to earth's atmosphere, ft/sec 
Vy vertical component of V, ft/sec 
y altitude, ft 

0 constant, l/22,000, l/ft 

T ratio of specific heat for air, 1.4 

€ time parameter 

Sg angle of flight path at reentrance with respect to horizontal, deg 
[x air viscosity, slugs/ft- sec 

p air density, slug/eu ft 

P 0 air density at sea level, 0.0034 slug/cu ft 
cp reentrant parameter, Oj-jA/ 20m sin 8g 

Subscripts : 

1 local 

max maximum value encountered along trajectory 
1 behind shock at stagnation point 


ANALYSIS 

The fluctuating surface pressures created by a turbulent bo undar y 
layer have been measured in a duct (ref. l) and on the surface of a wing 
(ref. 2) . In both cases the measurements were made by using microphones 
embedded in the surface and exposed to the airstream. Both sets of data 
showed that the root -mean- square pressure p at the surface varied 
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linearly with a dynamic pressure q of the flow. The ratio p/q was 
relatively independent of the Reynolds number or of the boundary-layer 
thickness. The data of reference 1 showed no effect of Mach number, 
whereas the flight data (ref. 2) indicated a Mach number effect only at 
low Mach numbers (below 0.55). Since the data of reference 2 were taken 
on an airfoil, the low Mach number effects could well result from angle- 
of-attack effects rather than from any real Mach number effect of itself. 
Unpublished data taken both in flight and in a large acoustic channel 
agree quite well with 1116 results of reference 1$ that is, no effect was 
noted of Mach number, Reynolds number, or boundary-layer_thickness on 
p/q. A study of all the available data indicates that p/q is approx- 

•z 

imately 4 to 5X10 for Mach numbers up to 0.8. For the case of a blunt- 
nosed reentry vehicle traveling at hypersonic speeds (M > 5) the local 
values of dynamic pressure may be quite large even though the local Mach 
numbers (behind the shock) are not high (less than 2.0). 

Since the value of p depends onJLocal quantities it would be ex- 
pected that the linear dependence of p on q would not be altered 
greatly by doubling the Mach number, that is, an increase from 0.8 to 1.6. 

Consider the case of a blunt-nosed body reentering the earth's at- 
mosphere at extremely high speeds as shown in the following sketch: 


At high Mach numbers an extremely strong bow wave precedes the body 
and the value of M^ is approximately 0.4 for a range of free-stream 
Mach numbers from 4 to 20. The flow behind the shock proceeding down- 
stream from the stagnation point speeds up and at some point on the body 
aft of the stagnation point the maximum value of the local dynamic pres- 
sure q^ outside the boundary layer is encountered. This maximum value 

in uniquely related to a particular value of M^ provided that the to- 
tal pressure of the flow around the body between the shock and the bound- 
ary layer remains constant and equal to the stagnation-point total 
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pressure P^- Then the maximum value of q^/P^ is 0.4312 and occurs at 
Mj = 1.415 (e.g. , ref. 3). It should be noted that the calculation pro- 
cedures used herein do not account for any effects associated with ex- 
tremely high temperatures such as dissociation or variable ratio of spe- 
cific heats. 


To calculate q.j at any point on the trajectory, it is only neces- 
sary to know the value of Pj_. In order to determine the maximum possi- 
ble value of P^_ for a given set of conditions, the following analysis 
was undertaken. The usual normal shock relations give 


£l 

p 



The maximum value of occurs at a particular combination of p and 

M, that is, at some particular point along the trajectory of the vehicle. 
Setting the derivative dP^/dy = 0 and solving for P-j_ will yield the 

maximum value of (P-, ) along the trajectory. The result of such a 

differentiation and algebraic rearrangement for Mq » 1 yields 


£ 

M dp 


-0.5 


Using the gas law p = pgRT and M = V/a, then, yields 


p dV 
V dp 


-0.5 


Allen (ref. 4) shows that the velocity at any point along the re- 
entrant trajectory can be expressed as 

V = Vge"^ 

where 

CpA 

^ 2£m sin Sg 

where Cjj is the body drag coefficient, A is the body area correspond- 
ing to Op, m is the mass of the body, and 8g is the reentrant angle 
of the body as It enters the atmosphere. The constant 0 appears in 
the atmospheric relation 
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The derivation of reference 4 assumes that C-q is essentially con- 
stant, that is, the parameter cp does not vary appreciably. Further- 
more, the values of 0g must be large enough so that gravitational ef- 
fects on the flight path are small . 


When these relations are used, the maximum value of P-^ occurs at 
pep = 0.5. Hence, for M >:> 1, 


(Pi) 


max 


0.1695 _2 
V-c 


cp r E 

and (q.7 ) mfty on the body corresponding to (P 1 ) max is given by 




max 


0.073 Vf 


The maximum q^ imm ediately behind the shock at the stagnation 
point is also obtained at the condition (Pi) and is given by 

, , 4 

(iiW - °- 01553 f 

The values of (q.j) max s-ud (q 1 ) m ax are unique functions of the al- 
titude y for specified values of cp (see appendix A) . The altitude y 
corresponding to (q^max or (l^max is P lc> ii et i as a function of q> in 
figure 1. It is evident that the peak dynamic pressures occur between 
40,000 and 95,000 feet for a range of 9 values from 1000 to 10,000. A 
plot of (q^ ) max as a function of the reentrant velocity for various val- 
ues of cp is given in. figure 2. 


The fact that the local dynamic pressures are high during reentry 
(500 lb/ft or greater) does not necessarily infer that a turbulent bound- 
ary layer exists on the vehicle. One important consideration that deter- 
mines the existence of a turbulent boundary layer is the Reynolds number 
associated with the flow. The relation for Reynolds number per foot 
(based on conditions immediately behind the shock at the nose 
= PpTp/ lj_) is shown in appendix B to be 



219X10 9 


max 


at (qj) 
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Furthermore , it can be readily shown (appendix B) that th e 
number per foot Re/ft based on free-stresm conditions at co ndi t 
corresponding to (q, ? ) TnRy is given by 

Re/ft = H X10 6 
<P 

Figure 3 shows a plot of Re-jy^ as a function of cp for reentrant 
velocities from 10,000 to 25,000 feet per Second. From the figure, it is 
evident that over the whole range of cp V E the value of Re^y^ 

varies by only a little more than a decade. Furthermore, even for a mod- 
erate size missile, values of Reynolds number well over one million w ill 
be achieved, which ind icates a high probability of the existence of a 
turbulent boundary layer. 

As pointed out previously, the value of (q? occurs at a local 

u max 

Mach number of about 1.4, and since p/q is essentially constant over a 
wide range of subsonic Mach numbers, it would be expected that there 
wouild be no large changes in p/q. at low supersonic speeds. There is, 
of course, the possibility that local shock-wave - boundary-layer inter- 
actions could greatly alter local values of p/q.. Such an effect is 
quite unlikely in the nose-cone region under discussion because the pres- 
sure gradients around the nose are favorable. 

It would appear, therefore, that assuming a constant value for p/q 

equal to 4. 5xl0 - ^ would_be justified for most preliminary calculations. 
Even with a change in p/q of 2, the sound pressure levels would vary 
by only 6 decibels. The sound pressures associated with (q, ) can be 

calculated using this assumption, and a graph of the ms-yirmim so und pres- 
sure levels SEL max as a function of 9 is given in figure 4 for var- 
ious values of initial reentrant velocity. Current values of cp «nri 
V E would give values of SFI^g^ of approximately 150 to 160 decibels 

for a long-range ballistic missile. It should be remembered that the 
values of sound pressure level will vary over the entire surface anil that 
SRLj^ax corresponds only to the maximum value that will occur only at one 

location aft of the stagnation point. However, since the magni t ude of 
the sound pressures is quite large compared with the us ual values it 
would appear that consideration should be given during the designing of 
long-range ballistics missiles to minimize the effects of these so und 
pressures properly. 
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One factor that would tend to alleviate both the external and inter- 
nal noise problems is the extremely short duration of exposure to the 
earth's atmosphere. The total time elapsed from the entrance of the nose 
cone into the earth's atmosphere until its impact on the surface would 
probably not exceed 120 seconds. Only during a portion of this time 
would the local dynamic pressures be high enough to result in high noise 
levels. Calculations of the variation of the maximum value of local dy- 
namic pressure with altitude for a range of cp values are shown in fig- 
ure 5. In each case the q^ value increases sharply below an altitude 

of 140,000 feet, reaches a peak value (o.j) max at some altitude between 

100,000 and 50,000 feet (dependent on cp) , and then decreases sharply as 
the altitude decreases. For values of cp between 10,000 and 3000, the 
high values of q^ exist between 140,000 and 50,000 feet. For lower cp 
values, the high values of q^ are encountered down to 30,000 feet. A 
method for calculating the time during which high values of q-j exist 
is given in appendix C. 

For the cases of current interest, the critical values of q j are 
encountered between 140,000 and 50,000 feet. The time for the missile 
to traverse this distance should give a reasonable estimate of the dura- 
tion of exposure. Figure 6 shows a time parameter s for the missile 
to fall from 140,000 to 50,000 feet as a function of parameter cp. The 
time parameter e is equal to ptVg sin 8g, where t is the actual time 

(sec) for the missile to traverse the given distance. For missiles of 
current interest the time is of the order of 25 seconds. For future mis- 
siles the times will probably approach approximately one-half this value. 


CONCLUDING REMARKS 

The results presented herein give an indication of the magnitude of 
the surface sound pressure levels that may be encountered during reentry. 
These levels are high (of the order of 150 db) but the time of exposure 
to the earth's atmosphere is short, that is, approximately 25 seconds. 
Furthermore, the Reynolds numbers are high enough to yield turbulent 
boundary layers. Current efforts to reduce heat transfer by maintaining 
large lengths of laminar boundary layers would greatly reduce any prob- 
lems that might be created by excessively high sound pressure levels. 

On the other hand, the use of ablation-type nose cones with their inevi- 
tably higher surface roughness should increase the order of magnitude of 
the sound pressures. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, May 8, 1958 
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APPENDIX A 

ALTITUDE AT WHICH MAXIMUM LOCAL DYNAMIC 
ERESSURE IS ENCOUNTERED 

The condition that yields Cq.2.)max or fa ls P9 = 0- 5 85 de ' 

rived in the ANALYSIS section. 

Hence, 

p 0 e" Py cp = 0.5 

where 

p = 1/22,000; p 0 = 0.0034 

Then 

e Py = 2p Q cp 
In e Py • In 2p Q cp 
Py = In 2p Q q> 

y = ^ In 2p 0 q> 

, = 22,000 ln(0.0068q>) 
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APPENDIX B 

REYNOLDS NUMBER PER FOOT BASED ON CONDITIONS IMMEDIATELY 
BEHIND SHOCK AT NOSE (p^/^) CORRESPONDING TO 
MAXIMUM LOCAL DYNAMIC PRESSURE 


Re 


l/ft 


P 1 V 1 pV 


= NT Slnce P 1 V 1 = pV 


h. M-i 


^ = 3. 44X10” 9 T°' 75 


t ‘h m 2 for M ’ 5 




* 3.44X10 


-9 [ JL (T\ rnl 
[_36 \ a/ J 


0.75 


Since the atmospheric relation p = p 0 e - ^ is an isothermal atmos- 


phere where T = constant = 392° R, 


*L “ 3 - 4 « L0 ' 9 [(?6)(^l) 2 392 ] ^ 


^1 


3 2.93X10' 


■12 ^1. 


Since pep = 0.5 for (<^)max> then V = V E e“ 0,5 = 0.606 V 3 


E 


Re 


l/ft 


pV 


5.415X10 11 P ^415)00 


11 


2. 93X10 " 12 V 1 * 5 


la). 


Vv 


VO. 606 V, 


E 


219X10- 

9VVg 


Furthermore, since V = 0.606 Vg at conditions for (<3.}) max and 
a = 971, then M = Vg/1600. 

Therefore, the free-stream Reynolds number per foot Re/ft is 
given by 

py ¥°- 606 7 k _ v e „ 

9 


Re/ft = R1 = 


V- 


3.03X10" 
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APPENDIX C 


CALCULATION OF TIME DURING WHICH HIGH VALUES OF 

LOCAL DYNAMIC IRESSURE ARE ENCOUNTERED 

In reference 4 it is shown that the vertical component of velocity 
Vy of a ballistic vehicle is given by 

V y = V sin % . $ 

-CDn 

where V = V-^e on all points of the trajectory. Hence, the t ime for 
a vehicle to cover the distance between two values of y is given by 


where 



t = 


-1 

PV E sin Qg 



dp 

pe^P 


t = 


1 

PV E sin Sg 




log p + p + 


Ccpp ) 2 

4 




Since p-j_ and Pg are unique functions of y, it is possible to calcu- 
late t for any chosen values of y-, and yg for particular values of 
<p, Vg, and sin 0g. If the time is expressed in terms of the dimension- 
less time parameter e = 0tVg sin then e is a function of cp only 
and can be shown by a single curve as indicated in figure 6. 
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Maximum local dynamic pressure during reentry 
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Figure 2. - Maximum local dynamic pressure around "body during reentrance as a function of re- 
entrant velocity. 
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7igure 3. - Reynolds number per foot corresponding to msud-inum local dynamic pressure as a function of re- 
entrant parameter . 
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Figure 6. - Time parameter for the missile to fall from 140,000 
to 50,000 feet as a function of reentrant parameter. 


NACA - La.igley Field, V*. 



